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Abstract. In this topical review we discuss the connections between chaos, decoherence 
and quantum cosmology. We understand chaos as classical chaos in systems with a 
finite number of degrees of freedom, decoherence as environment induced decoherence and 
quantum cosmology as the theory of the Wheeler - DeWitt equation or else the consistent 
history formulation thereof, first in mini super spaces and later through its extension to midi 
super spaces. The overall conclusion is that consideration of decoherence is necessary (and 
probably sufficient) to sustain an interpretation of quantum cosmology based on the Wave 
function of the Universe adopting a Wentzel - Kramers - Brillouin form for large Universes, 
but a definitive account of the semiclassical transition in classically chaotic cosmological 
models is not available in the literature yet. 
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1. Introduction 

The goal of this topical review is to explore the connections between the concepts of chaos, 
decoherence and quantum cosmology. Each of these has a long history and development 
independently of the other two, and the relationships between chaos and decoherence are 
varied and deep. When it comes to their relationship to quantum cosmology, on the other 
hand, it is fair to say that the available literature barely scratches the surface of the subject. 
Therefore our goal is not only to review what has been done to date, but also to show areas 
which remain underdeveloped, sometimes in spite of the fact that the necessary tools are 
already available in non-quantum cosmological contexts. 

Two warnings to the reader are called for. First, we shall make no attempt to 
be systematic in our presentation of chaos (TOj [161 E21 ED 1248] 1256] and decoherence 
|1211 12531 1275] , we shall only pinpoint those aspects of these subjects which have been found 
to be relevant to quantum cosmology. Most importantly, we shall neither attempt a complete 
presentation of quantum cosmology itself [35j [36j I182|. I2T0|. 1271] . For the purposes of this 
review, quantum cosmology means the theory of the Wheeler-DeWitt equation restricted to 
Friedmann - Robertson - Walker mini super spaces (to be defined below), later enlarged to 
midi super spaces by including also linearized primordial fluctuations (whereby the former 
mini super spaces come to be regarded as quantum open systems). We have made this choice 
not because we think this is the most important or promising brand of quantum cosmology 
(if asked point blank, we would rather give that award to the loop quantization program 
p^2| [13] ) but because that is the framework in which the connections of quantum cosmology 
to chaos and decoherence have been most systematically addressed, and where the language 
remains closest to the other area where cosmology and quantum physics intertwine, namely, 
the generation of primordial fluctuations during inflation [731 1225] 1226] . 

Let us try to put the above choice of subject matter in perspective. 

It is certainly not coincidence that after a strong paradigm of quantum field theory 
based on the renormalization program emerged in the fifties and early sixties, people started 
to worry about extending it to other areas of physics, such as condensed matter theory 
[90] EH [mj EM], turbulence [THTJ [188] and gravity [Ml QUE]. In the case of gravity, which 
is a constrained system development had to await for the unraveling of the Hamiltonian 
structure of the theory [9?]. Once this was accomplished by Arnowitt, Deser and Misner 
(ADM) [HI 1217] . the road was open for a systematic quantization of gravity on the lines of 
Dirac's quantization [951 ESI EH 1138] . This program was undertaken by Wheeler, DeWitt 
and collaborators, and resulted in the formulation of the Wheeler - DeWitt equation, which 
is the cornerstone of our subject [921 [931 12761 [1901 HSU- 

This first wave of quantum cosmology research stumbled upon the problem of actually 
finding solutions to the constraints of the theory (namely, states). Rather than making no 
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progress at all, Misner and others hit upon the idea of working on (very) restricted models 
where the momentum constraints where trivial, and only the Hamiltonian constraint, or 
Wheeler - DeWitt equation, remained as dynamical law [271 EQj ELI EE31 EH EES EEEJ E2] . 
At first, these so-called mini super space models were offered as just mathematically 
consistent models which shared one important property (time-reparametrization invariance) 
with full quantum cosmology, and could therefore be of use to explore outstanding issues 
such as the emergence of time and the correct interpretation of the theory. Eventually mini 
super spaces came to be regarded as all that and also as some kind of Born - Oppenheimer 
approximation to the full theory, so that the study of minisuperspaces was directly relevant 
to cosmology |192l [193J. The problem of actually solving the constraints of the theory was 
reopened much later by the loop quantization program, which we shall not discuss here 

pi Eg. 

At about this point the program developed by Wheeler, Misner and their associates 
merged with a powerful current coming from the (then) Soviet Union, where Landau, 
Belinskii, Lifshitz, Khalatnikov and their collaborators had made foundational work on 
the fluctuations of Friedmann - Robertson - Walker cosmologies, and described the 
homogeneous Bianchi-type cosmologies, including the notorious Bianchi type IX model 
(2HE51I251 0231 023 Q2HI 0231 USE]. In time, Bianchi IX came to play a role, with respect 
to cosmology and chaos, not unlike the one played by the Schwarzschild Black Hole with 
respect to gravity and thermodynamics. 

A third moment in the development of quantum cosmology came when the emphasis 
shifted from the analysis of the Wheeler - DeWitt equation in general, to the identification 
of which particular solution to that equation described the state of our Universe, and the 
correct way to derive cosmo logical predictions from it. Here the foundational work is no 
doubt that of Hartle and Hawking, which brought the formalism of euclidean path integrals 
to bear on the subject [1331 EE1 EM EM EM EM ED], and of Vilenkin [2721 [273], who 
proposed a boundary condition for the wave function of the Universe allowing a direct link 
to inflationary cosmology. 

If at the time of Misner's Mixmaster cosmology the main issue was how our Friedmann 
- Robertson - Walker model could have emerged from such chaotic, anisotropic early stages 
[221 E31 EH [2081 E091 EES EEE 1213 U, after the inflationary paradigm became dominant 
that question (if not answered [511 E31 [761 27]) l° s t poignancy, since it was accepted that as 
long as there were one single inflationary patch in the whole Universe that could account for 
our existence. Thereby the relevant question became not only the properties of the classical 
Universe emerging from the primordial quantum state of the Universe, but also how such a 
classical Universe could have emerged at all. 

This was the first point of contact between chaos and quantum cosmology, because, as 
in quantum mechanics more generally, a classical Universe was identified with a Wentzel 
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- Kramers - Brillouin (WKB) solution to the Wheeler-Dewitt equation [133J, and it is 
known since Einstein that nonintegrability set strong limits to the applicability of the WKB 
approximation [104] . It soon became clear that if classicality was identified with WKB form, 
then an Universe could be classical near the cosmological singularity and turn quantum upon 
expansion, quite the opposite to expectations [6H [891 H95L 1233] . 

In ordinary quantum mechanics, it is generally accepted that properties of a system 
become actual when they are measured and recorded. This presupposes there are measuring 
and recording devices external to the system itself, and whose behavior is classical in at least 
some important respects. Clearly such a framework is inappropriate for cosmology, where 
the system encompasses the whole Universe. This led Hartle and Gell-Mann to develop from 
scratch a quantum theory of closed systems, by adopting the consistent histories approach 
to quantum mechanics [Ml [1021 [1031 EH EH EH [1201 EU E21 E21 Ell E3 [1531 ESI 
I176[ 1231] . This meant to concentrate not on the wave functions of the Universe at a given 
point in cosmological evolution, but rather on full histories of the Universe and their mutual 
consistency. 

Same as in quantum mechanics at large, the emergence of a classical state from 
the quantum maelstrom involves a process of loss of quantum coherence, the so-called 
decoherence. Since we will take as axiomatic that the quantum mechanics of a closed system 
is unitary, if decoherence occurs it means the system is not closed. It is interacting with an 
environment, and by simplicity we shall consider only cases where the environment is actually 
a physical system on its own right. More abstract settings, for example the environment being 
high order correlations of the system itself, are conceivable [H [671 ES]; we shall not consider 
them here for pedagogical (and space) reasons. This leaves us squarely within the so-called 
environment induced decoherence paradigm [SHI El EH [2661 12581 (2821 [2831 1283 [2851 12861 
1287] . 

Actually, environment induced decoherence explains how an ensemble of classical 
Universes emerges from the quantum state. There is a second part to the question, which is 
how a single classical Universe is then selected from the ensemble. This process has many 
points in common with spontaneous symmetry breaking and also with quantum measurement 
theory, and we regret we have nothing new to offer beyond what is already known in those 
fields pQ. 

The goal of the environment induced decoherence program is not to explain why the 
density matrix of the system (in our case, the system is the Universe, or rather its projection 
upon the mini super space) becomes diagonal, since the density matrix is always diagonal 
in some basis, but rather to point out some preferred basis, and to show that the density 
matrix inexorably becomes diagonal in that basis. 

The problem of choosing one single preferred basis out of the bewildering array of 
possibilities (in the consistent histories program there is the equivalent problem of which 
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variables should be monitored as part of the description of the consistent history set) is not 
particularly new in Science. Max Weber confronted essentially the same problem in trying 
to find a consistent way of describing capitalism |274] . Weber's answer was that the selected 
basis should answer to certain internal consistency requirements (in our case, for example, 
the predictability sieve |245[ 1288] to be introduced below provides such a consistency check), 
but should reflect our interests also - the choice of preferred basis is informed by exactly 
what we are trying to describe. Weber's criteria are also good for quantum cosmology. 

Once we accept that our problem is the emergence of a classical ensemble of 
Universes with respect to a restricted set of properties chosen a priori, then environment 
induced decoherence has much to say about the conflict between the semiclassical limit 
and nonintegrability we have found above. Indeed, not only the validity of the WKB 
approximation is restored, but also the process of decoherence is largely universal, that 
is, environment independent |288[ 1289} I291[ [290] . We shall give a precise description of what 
we mean below. 

Surprisingly, however, the consequences for quantum cosmology of these insights from 
environment induced decoherence have not been worked out yet to our knowledge. This will 
be the point where we part with the reader, pointing out the new lands we could quickly 
conquer, if we set to work on it. 

This paper is organized as follows. In Section 2 we provide the necessary background 
on classical and quantum Hamiltonian cosmology, and the consistent histories formulation 
of the latter. In Section 3 we discuss classical cosmological chaos, both in Bianchi IX 
and in Friedmann - Robertson - Walker models, and how the fact of classical chaos raises 
some foundational issues for quantum cosmology. In Section 4 we discuss decoherence in 
cosmology, by regarding the mini super spaces considered so far as quantum open systems 
embedded into an environment of gauge invariant primordial fluctuations. We conclude 
pointing out a few highlights from the literature. 

2. Quantum Cosmology 

2.1. Hamiltonian Cosmology 

The first step in implementing the quantum cosmology program is to formulate General 
Relativity in a language adequate for its quantization. This language is the canonical 
formulation pH EGl l25T|l252]. 

We begin by writing the so/called ADM decomposition of the interval element 

ds 2 = gffldx^dx" = -N 2 dt 2 + g® [dx a + N a dt) (dx b + N b dt) (1) 

We follow Misner-Thorne- Wheeler (MTW) conventions throughout |217j . indexes /i, v run 
from to 3, indexes a, b run from 1 to 3 and x° = t. We also assume natural units with 



Chaos, decoherence and quantum cosmology 



h — c — 1. g^ b is the induced metric on t = constant surfaces. 
The Einstein - Hilbert action reads 

S 



(2) 

plus eventually a cosmo logical constant, m p is Planck's mass and we have discarded a total 
divergence. Further integrations by parts allow us to reduce this to 



S 



d"x NJg 



K ab K ab 
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Of 



(3) 



(4) 



is the extrinsic curvature of the t = constant surfaces and K = g^ b K ab . The symbol 
denotes a covariant derivative with respect to g 



(3) 



ab- 



The action |2] does not depend on the time derivatives of the lapse N or shift N a . These 
are Lagrange multipliers enforcing the constraints of the theory. The canonical momenta 
conjugated to g^\ b are 
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The Hamiltonian density is (discarding a total divergence) 
H = NH° - 2N n H a 



where 



1 



m*\/g^ 
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rnt\/ g 
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(5) 
(6) 
(7) 



7T 



(3>6 



-"" = g^ ab 7T^ b and is the intrinsic scalar curvature computed from g^\ b . The 



-(3) 



constraints of the theory are H° = H a = 0; the H° and H a are the generators of coordinate 
transformations. HH° is the so-called Hamiltonian constraint, while the H a are the momenta 
constraints. 

The space of all three metrics is the so-called super space [9U 1111] . The kinetic term 
in the Hamiltonian H° may be written as proportional to (1/2) G^^n^ahTc^cd with the 
super space "metric" 



Q(ab)(cd) 



(3>c (3>d _j_ (3)od (3>c 



g (3)ab g (3)cd 



(9) 



Point by point, this metric has signature 



+ + ++), with the "time-like" direction 



corresponding to the conformal degree of freedom of the three metric, as we shall show 
presently 
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2.2. Singling out the conformal degree of freedom 

A model where matter is represented by a single conformally coupled scalar field is possibly 
the simplest one which nevertheless is rich enough to allow for the discussion of chaos and 
decoherence. It will therefore be our model of choice to illustrate the different points to 
follow. With this in mind, we begin by formulating general relativity in a way where the 
dynamics of the conformal degree of freedom is singled out. 

Let us write the four dimensional metric as gj$ = a 2 g^.The three dimensional metric 
shall be g^\b = a 2 g a b- a = e n is a generic scalar field, admits an ADM decomposition 
as above, with lapse N = aN and shift N a = N a . To make this decomposition unique we 
should enforce some condition on g ab , such as constraining the determinant g constant. After 
splitting variables in the 3 + 1 way and separating out the conformal degree of freedom we 
get the connection as 



F (4> = 
fit/ 

The Riemann tensor 



(10) 



The Ricci tensor 

R (4 lu = Rau - g,ug TX n. M - (d - 2) Q, av + {d-2) 

and the scalar curvature 

R^ = 1 |i? - (d - 1) (2V 2 VL + (d-2) (Vfi) 2 ) | 

where 



(11) 



(12) 



(13) 



v 2 n = 



V ^ [g^ - n"n"] Ny/gQ v 



N^g 

(Vttf = [g^ - n^n"] ^, „ 



(14) 



Of course, here d = 4; setting d = 3 we obtain the formulae appropriate to the intrinsic 
curvature in the t = constant surfaces. 
The extrinsic curvature 
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The action 



K ab = a ^K ab - -1 [Q ta N b + tt b N a + g ab (ft - iV e ft e )]j 
= aK ab - [a a iV 6 + a b iV a + g ab (a - N e a, e )] 



(15) 



S = m 2 , j rf 4 a; iV^a 2 ji? - ^ 2 + X bo X ba - A (n t _ 7V a ft a - ^ 
+ m 2 / rf 4 x 2iY^ {- Aa 2 + 3<7 a6 a (16) 



The canonical momentum conjugated to a is 



12m 



N 



n, - N a n n - 



NK 



While the canonical momenta conjugated to the g ab are (recall that g ab g a b,o = 0) 



n ab = -a'mt^g- 



2_2 



K 



ab 



-g ab K 
3 y 



(17) 



(18) 



As expected, the kinetic term associated to the conformal degree of freedom appears 
with a negative sign. 

2.3. Adding conformal matter 

We now add to the above model a conformal scalar field $. The action is 



s m = 4l d * x A + [™ 2 + ^TT) R{A) \ $ 



(19) 



where d — 4. We write $ = a ( d 2 )/ 2 $ and proceed to make explicit the conformal degree 
of freedom as above to get 



$2 S 
12 



1 

3" 



R-^K 2 + K ba K ba 



S= J d A xN^g-(mla 2 - 

- 6m 2 , J d A x N- 1 ^ (a, t - N a a >a - l -NKa 
+ 2m 2 , J d A x N^{-Aa 2 + 3g a \ a a, b } 

+ lf d 4 x N~ l ^g- t - jV°$ a - \nK*} 
d A x Ny/g 



(20) 
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So now the canonical momenta conjugated to g a b read 
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The canonical momentum conjugated to a remains the same 

NK~ 



12m 2 p r 



n, t - N a n, a 



and the field canonical momentum is 



p = N- 1 ^ ( $, t - N a $, a - -NK$) 



leading to the Hamiltonian density 
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2.4- The Wheeler- DeWitt equation 



(21) 



(22) 



(23) 



(24) 



We proceed now to quantize the theory along the lines of Dirac's quantization for constrained 
systems [951 1961 1981 [138] . 

Leaving aside the subtleties associated with the super space geometry, the basic idea is 
to associate the super space "coordinates" g ab by the operators "multiplication by g a b" , and 
the momenta 7i ab by the variational derivative with respect to g a b- After sorting out operator 
ordering ambiguities the constraints are mapped into operators. The states of the theory are 
then the joint null eigenstates of the constraints. Particularly, the Hamiltonian constraint 



H n V = 



(25) 



takes the form of a time-independent Schrodinger equation for the so-called "wave function 
of the Universe" This is the Wheeler - DeWitt equation [921 [931 1223 [Mil EE]. 
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It is natural to seek an operator ordering such that the Wheeler - DeWitt equation ends 
up being covariant in super space, if it is understood that $ is itself a scalar. If we introduce 
indexes A, B running from to 5 to represent all independent pairs of indexes (ab) and write 
x A = g a b then the Wheeler - DeWitt equation becomes a wave equation is super space 



where G AB = G {ab){cd) is the super space metric, G < its determinant and R its curvature. 
The undetermined constant £ is a remnant of operator ordering ambiguity. In the following 
we shall assume "minimal" coupling £ = (see however |131[ 1214] ) . 

Actually, it is very hard to solve the constraints of the theory in these variables. This 
has led Ashtekhar and others to take on the problem from a different viewpoint, eventually 
leading to the loop quantum cosmology program. We shall not discuss these developments 
here. 

2.5. The Hartle- Hawking and Vilenkin boundary conditions 

Given that the Wheeler - DeWitt equation is analogous to a Schrodinger equation is super 
space, Hartle and Hawking have proposed to solve it by importing the euclidean path integral 
methods of ordinary quantum theory. Their ansatz reads |133[ 11371 1146[ 1147] 11491 H50j 1151] 



Here g a b is a three metric defined on some Cauchy surface S. The integral is over all four 
dimensional non- singular Riemannian manyfolds having E as their boundary, and over all 
euclidean metrics on them such that g a b is the induced metric on E. Observe that while the 
physical space is suppossed to originate from a Big bang, which therefore defines a boundary 
to the past, the euclidean metrics in the Hartle - Hawking proposal have no boundary other 
than E. Se is the Einstein - Hilbert action for the euclidean metric. 

The integral in [27] is ill defined because Se is not nonnegative definite. This is related to 
the signature of the super space metric G AB . One possible way out is to perform a "second" 
Wick rotation (the first one was going from Lorentzian to Riemannian metrics), that is, 
one integrates over eucliean metrics with complex conformal factors. \I/ [g a b\ then is defined 
through analytical continuation, and may become complex in the process [133J . 

Vilenkin has proposed an alternative to the Hartle - Hawking where the Universe is 
created from nothing through a quantum tunneling process [2721 1273] . 

2. 6. Mini Super Space models 

Mini super space models are Hamiltonian systems that result from parametrising the metric 
and matter fields in the Universe leaving only a finite number of time-dependent but space- 




(26) 




(27) 
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homogeneous degrees of freedom indetermined. The parametrisation should be such that 
introducing the restricted metric in the Einstein - Hilbert action and then taking variations 
with respect to the parameters should be equivalent to introducing the restricted metric in 
the Einstein equations themselves. Observe that spatial homogeneity does not mean that 
the fields are independent of the space coordinates, but only that there is a group of motions 
which acts transitively on space and leaves the restricted metric invariant |214j . 

In a typical mini super space model the parametrization is chosen in such a way that 
the momenta constraint hold identically at the classical level. If the q a are the independent 
parameters, then the Einstein - Hilbert action becomes 



S = m z p I dtN 



2iV 2 ' 



(28) 



We have the freedom to choose N, which could depend on the q a . q denotes a time derivative, 
of course, only after we have chosen N we will know what kind of time is t. 
The canonical momenta are introduced in the usual way 

2 

Pa = -j^fapf (29) 

leading to the Hamiltonian constraint 

#o = + m 2 p U [q] = (30) 

The theory is quantized by replacing q a by the operator q a - and p a by (— i) d/dq a . Operator 
ordering is chosen in such a way that —f a/3 p a p/3 becomes the D'Alambertian associated to 
the metric f a /3- 

It is clear that simply ignoring the degrees of freedom other than the q a is, at least, 
a violation of Heisenberg's principle. However, if it can be shown that the q a are in some 
way the "slow" degrees of freedom, then the mini super space may be regarded as the 
Born - Oppenheimer approximation to quantum cosmology. Eventually, the theory may be 
improved by including the zero point energy of all other degrees of freedom into U [q] . We 
shall return to this below. 



2.7. The interpretation of QC 

We shall now discuss different proposals for an interpretation of the Wave Function of the 
Universe EH H32J [1391 UM US IT821I273] . 

Of course, the simplest interpretation is just to assume Born's rule: let |^| be the 
probability density for three geometries in super space. One obvious problem is that \I/ is 
not usually square integrable; however there are ways to factor out this divergence and come 
up with a viable normalization [136J. 
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A more conceptual problem is the following. We are used and intend to think of 
cosmology in terms of evolution. We wish statements such as "the Universe becomes classical 
when such a such occurs" to make sense. The Born rule would allow us to say than a smooth 
isotropic Universe is more or less likely than a foamy fractal one, but gives us no clue 
about which comes first. Of course there is no "time" strictu sensu in quantum cosmology, 
but the fact itself of the Lorentzian nature of the super space metric suggests there is a 
causal structure, with well defined Cauchy surfaces arranged in increasing order of some 
"Heraclitean" variable |267j . Indeed, for most practical purposes, the volume of a spatial 
three section, or the average scale factor for open Universes, is a good enough ersatz for 
"time" . 

The analogy of the Wheeler - DeWitt equation with the Klein - Gordon equation 
moreover teaches us how to construct a current whose flux over Cauchy surfaces is conserved. 
This suggests to normalize \I/ to unit total flux, and to consider the normal component of the 
current as a measure of relative probability within a given Cauchy surface. For the De Sitter 
mini super space above this is just the ordinary probability current (— i) (\]/*\l/ )a — ty^* a ). 



As we can see, it is not nonnegative definite, and therefore the probability interpretation is 
untenable. Even if we restrict ourselves to "positive frequency waves" for which the normal 
component of the current is nonnegative, the set of allowed wave functions will be different 
for different Cauchy surfaces, a very unsatisfactory state of affairs. 

The drawbacks of these (and other, more sophisticated ones) interpretations of \I/ suggest 
than something drastic may be called for. Basically, we may take as an observational fact 
that the Universe became classical as it expanded, and in many examples \l/ adopts a WKB 
form for large values of the scale factor. If we merge these two insights, we obtain the 
"WKB interpretation" of quantum cosmology, namely, a wave function of the WKB form 
means that an observer within that Universe will perceive around her a classical Universe 
evolving under the particular solution of the Einstein equations such that p = S >a . [133] 

The situation when there is a coherent superposition of different WKB solutions is much 
more unclear [T8T] . It seems some degree of decoherence is necessary to turn such a pure 
state into an ensemble of WKB alternatives, each of which can then be interpreted of the 
prediction of a particular classical Universe. The square absolute values of the coefficients in 
the original superposition give the relative probability of each Universe within the ensemble. 

Most importantly, and finally getting squarely into our subject matter, the WKB 
interpretation will come under fire when the underlying Hamiltonian system is not integrable, 
because then \1/ may not take a WKB form at all, certainly not for large Universes 
[611 [891 11951 1233] . It is here that we shall see that only a joint consideration of chaos 
and decoherence yields a viable quantum cosmology 
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2. 8. Hartle Gell-Mann formalism 

To analyze this question we shall adopt the consistent histories approach to quantum 
mechanics, in the version advanced by Gell-Mann and Hartle [lUlj 1102} 11031 H171 11181 ESI 
[1201 EH E3 QUil ES ESI [1531 EE! ESI ED]. The idea is to define a history by a set 
of projectors P a acting at times tj. In canonical terms, a history is given by an evolution 
of the state vector such that at every time t iy it belongs to the proper space of P a (tj). In 
path integral terms, the projectors are represented by window functions w a [x (£«)], which 
take on unit value if the instantaneous configuration x satisfies the requirements of the 
history a, and vanish otherwise. The limiting case of a fine-grained history, namely, when 
x (t) is specified for all times, is assigned an amplitude expiS/h, as usual in the Feynman 
path integral formulation. The amplitude for a coarse-grained history defined by window 
functions w a [x (ti)} is defined by the superposition 

A[a] = jDxe^[x(0)]^llw a [x(U)]^ (31) 
The probability is naturally expressed in terms of a closed time path integral 

P[a] = J DxDx' e^-^/Vtx^y^ln^tx^lUn^a^ (32) 

In this way we may assign a probability to any coarse-grained history, but these probability 
assignments are not generally consistent, namely, the probabilities of two mutually exclusive 
histories do not generally add up. Indeed, let us define the decoherence functional of two 
histories a and j3 

V [a, 0\ = j DxDx' jV-Wp [x (0) , x' (0)] ^[w a [x (U)] j U[ w p [x' fa)) J (33) 

P [a] = V [a, a] but P [a V 0\ = V [a, a) + V [(3, (3} + 2Re£> [a, (3} ^ P [a] + P [/?]. The 
probability sum rule P [a V 0\ = P [a] + P only applies when the third term vanishes, 
and in particular when there is strong decoherence, V [a, (3} = for a ^ (3. As physicists, who 
deal with reality, we shall be satisfied that a set of mutually exclusive histories is consistent 
when \T> [a, f3] \ <C V [a, a] , V [(3, (3] whenever a ^ (3. 

To a certain extent, the consistent histores approach has the WKB interpretation built 
in, because the probability of any history will be negligible, at least in the limit m p — > oo 
unless the history contains stationary points of the Einstein-Hilbert action. 

A simple set of consistent histories refers to the values of conserved quantities |148j . 
First observe that the path integral expression fl33|) translates into the canonical expression 



V [a, (3}=Tr If 



T 



p(0) . (34) 
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The projectors at different times are related in the usual way P a (t) = U (t) P a (0) W (t) . 
If a projector commutes with the Hamiltonian, then it is time- independent, and expression 
( |34|) collapses unless all projectors are indeed identical. The only histories with nonzero 
probabilities are those defined by ranges of conserved quantities in the initial state, and they 
are automatically consistent if these ranges do not overlap. 

3. Chaos and Cosmology 

The paradigmatic chaotic dynamical system is probably the Baker transformation [HH] . This 
is a map of the unit square into itself. If (x, y) is a point in the square, then it is mapped 
into (2x — [2x] , (y + [2s]) /2), where [] stands for integer part. Alternatively, if 

oo 

x= Ec^— 1 

i=0 
oo 

y=Y. c * 2 ~ l (35) 

i=l 

with every q = or 1. Then the original point is represented by the doubly infinite sequence 
Q = q , and it is mapped into the sequence q = cf\. The mapping is obviously time- 
reversal invariant and area preserving. 

There are two typical ways to quantify the chaoticity of the Baker transformation. 
Suppose we have two points a distance dx apart. Then each iteration will increase this 
distance by a factor e A , where A = In 2. We say the map has a positive Lyapunov exponent. 
Since it is area preserving, it must have a negative Lyapunov exponent too, in this case, for 
points separated along the y direction. 

Now suppose we measure the Cg components over a large number of iterations. No 
matter how many times we do it, we shall always be one bit of information short of predicting 
the next outcome. We say the map has a Kolmogorov entropy Sk — In 2. We notice the 
Kolmogorov entropy is just the same as the (only) positive Lyapunov exponent. This is a 
particular instance of the so-called Pesin theorem. 

The Baker transformation is discontinuous at the x = 1/2 line, but this is nonessential. 
We could as well stretch the square and fold it back into itself without cutting it, obtaining 
the so-called Smale horseshoe. This way many points in the square are not mapped back 
into it, but there is an uncountable set of points that belong to the square and remain in it 
after any number of iterations, either back o forth. Restricted to this set, the dynamics is 
quite the same as the original Baker. 

The amazing thing is that the phase space of nonintegrable Hamiltonian systems 
generally contains subsets where the dynamics reproduces the Smale horseshoe [TQl I197[ [246J . 
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3.1. Chaotic cosmological models 

As we said in the Introduction, mini super spaces are usually restricted to homogeneous space 
times. That still leaves open the alternative that space may be isotropic (and thus the model 
belongs to the Friedmann - Robertson - Walker (FRW) class) or anisotropic (thus falling 
under one of the nine types of the Bianchi classification). We shall be mostly concerned with 
the first (FRW) CctSG, clS it is the most directly relevant to inflationary cosmology. However, 
nontrivial dynamical behavior was first found in the Bianchi IX type mini super space, 
and this beautiful problem still is a major motivation in the development of the subject. 
Therefore, we shall point out some few entry points to the literature on Bianchi IX models 
before resumming our analysis of FRW ones. 

3.1.1. The Bianchi IX models The metric of a homogeneous space time may be written as 
ds 2 = -N 2 (t) dt 2 + g ab (t) a a <g) a b (36) 

The cx a 's are three independent one-forms obeying 

da a = C a hc a h A o c (37) 

Homogeneity shows in the C£ c being space-independent. The C£ c are antisymmetric in (b, c) 
but may be of mixed symmetry with respect to a. Classifying them according to their 
symmetry we obtain the nine Bianchi types. The most interesting for our purposes are 
Bianchi I, where all C£, = 0, and Bianchi IX, where C£ c = e a bc is totally antisymmetric. 

In the Bianchi I case, also known as Kasner Universe, we may simply take a a = dx a . 
Then the solution of the vacuum Einstein equations reads 

ds * = _ R 2 {t) dt 2 + t 2 Ql ^iy + t 2 Q2 ^ 2 + t 2 Q3 ^ 2 (3g) 

The exponents satisfy qi + + Q3 — Q 2 + Q2 + Q3 = 1- They may be parameterised in terms 
of a single parameter u ranging from 1 to 00 as 



Q2 



— u 

1 + U + V? 

1 +u 

1 + U + V? 

u + u 2 



* = n 7 2 ( 39 ) 

1 + u + u z 
or some permutation thereof. 

The Bianchi IX evolution may be understood as a sequence of "Kasner epochs", 
themselves grouped within "Kasner eras". In each Kasner epoch, the metric may be 
approximated by [2H1 (with the appropriate one forms, of course) and therefore characterized 
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by a parameter u with some well defined value. If u > 2, then the Kasner epoch eventually 
evolves into another Kasner epoch with u' — u — 1. However, if 1 < u < 2, then the Kasner 
era ends, and another era begins with u' = 1/ (u — 1). 

It follows that if u n is the value of u at the first epoch of the n-th era, then the number 
of epochs in the era is [u] — 1 ([] standing for integer part). Let us "compactify" the w n 's 
by writing x n = u~ l . Then < x n < 1, and each x n is retrieved from the preceding one 
through the mapping 



This mapping has a positive Kolmogorov entropy, and thus it is strongly chaotic. To the best 
of our knowledge, this insight was the first quantitative indication of chaos in cosmological 
models [261 1201 1213] . 

The Bianchi IX Universe has been subsequently analysed by a large variety of techniques, 
confirming the original insight [191 EU E3J EHl EH1 EH SHI S31 SSI SHI [85l EZl [HHl EEl [EE3l fTTll 
fT52l [TSil H55l 12181 12161 12291 12TTI 12^91 [2501 12631 [280] . 

We are not aware of a direct proof of Bianchi IX chaos in terms of Ashtekhar variables. 
Bianchi IX may be regarded as a perturbation of Bianchi VIII, which is explicitly integrable. 
However, when one tries to carry the computation from Bianchi VIII to Bianchi IX, 
perturbation theory breaks down. This is of course in agreement with expectations, but 
falls short of an actual proof of non integrability [123[ 178] . 



%n— 1 



3.1.2. The Friedmann- Roberts on- Walker models A Hamiltonian system with n degrees of 
freedom is integrable if it has n constants of motion in involution, which means that all 
the Poisson brackets between any two of these constants vanish. A FRW mini super space 
has only one geometrical degree of freedom, the radius of the Universe a, and at least one 
constant of motion, the Hamiltonian itself. Therefore, to be non integrable it must contain 
non gravitational degrees of freedom, i. e. matter fields, as well. The simplest FRW mini 
super space model where chaos may be investigated has, besides a, a single homogeneous 
scalar field (f). It is extremely suggestive that already this basic setup leads to nontrivial 
dynamical behavior [371 E2J [58] . 

The model investigated by Calzetta and El Hasi is spatially closed; the matter field is 
free, massive and conformally coupled. The model describes a periodic Universe undergoing 
a sequence of big bangs and collapses. Non integrability becomes apparent when the whole 
sequence is considered |224j . rather than a single cosmic cycle. More realistic models 
including more matter fields already show observable consequences of chaos within a single 
cosmic cycle [591 EH] . See also |230[ 1265] . 

By now the consensus is that FRW Universes filled with an scalar field are non integrable 
except for a few exceptional cases. The field may be real or complex, and conformally or 
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minimally coupled, and there may be a cosmological constant or not [531 12011 [202J 

Let us discuss briefly the model analyzed by Calzetta and El Hasi [58J. It assumes a 
FRW spatially closed geometry, that is, a line element 

ds 2 = a 2 { V ) [-d V 2 + d X 2 + sin 2 X {d6 2 + sin 2 6dtp 2 )} (41) 

Where < <p < 2ir, < 8 < n, < x < 7r, and rj stands for "conformal" time. For 
concreteness, we shall consider only models starting from a cosmic singularity, that is, we 
restrict r\ to be positive, with a(0) = 0. Also, as indicated by the dynamics, we shall assume 
that after the Big Crunch ( that is, when a returns to ), a new cosmological cycle begins, 
now with a < 0. Therefore, a complete periodic orbit describes the birth and death of two 
Universes. The scalar curvature R = 6 and the conformal volume of an spatial surface is 
v = 2it 2 . It is convenient to rescale the conformal factor a = a s /y/2A7im p , the homogeneous 
matter field <p — 'Ps/ V2tt and the field mass m = \f2Amm v ii. Calling P s and p s the canonical 
momenta conjugated to a s and <fr s , respectively, we get the Hamiltonian 

H = {\)[-{P 2 s + a]) + (p 2 s + + /i 2 a 2 2 ] (42) 

The Hamiltonian constraint is H = 0. We introduce the "adiabatic" amplitude and phase j 
and "d 

(43) 




p s = J 2uj cos ■d (44) 

Where u 2 = 1 + fi 2 a 2 is the instantaneous frequency of the field. This transformation is 
canonical if we change the geometrical momentum from P s to P±, according to 

p ' = Pl + 2(iT^!) sm2tf (45 > 

The Hamiltonian becomes H = —(H + 5H), where the unperturbed Hamiltonian 

Ho = {\)[Pl + a 2 s }-ju (46) 



and the perturbation 



[i 2 a s P X 3 . 
on = — sin 2v + 

2cu 



A< a s J 
Au 2 



cos 4tf) (47) 



The dynamics of a s , as generated by H , is obviously bounded. The point a s = is a 
fixed point; it is stable if fi 2 j < 1, , and unstable otherwise. In this second case, there is 
an homoclinic loop associated with it. However, this orbit does not satisfy the Hamiltonian 
constraint; rather, we have Ho = —j on the homoclinic loop. 
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The equations of motion are simpler if written in terms of a new variable X — u, rather 
than a s itself. The transformation is canonical if we associate to X the momentum 

In the large j limit, with H ~ and ji fixed. X ^> 1 for most of an orbit, and the 
Hamiltonian simplifies to 

Ho = (^)Pl + {±)[(X - /i 2 j) 2 - (/i 2 j) 2 - 1] (49) 
If we parametrize 



P x = - — cos 2a (50) 
X = n 2 j + n\/~K sin 2a (51) 



H =[]-)[K- { M f - h (52) 



then 

We have succeeded in integrating the unperturbed motion, only that, because the last 
reparametrization involves j, the angle canonically conjugated to it is no longer but 

6 = tf- ii\[~K cos 2a (53) 

The perturbation causes resonances between the variables a and 6 and thus destroys the 
integrability of the system. Resonances occur whenever ojj ± tiuk = 0, io^k being the 
frequencies associated to 9 and a, respectively. We have Uj = fi 2 j and ujk = —1/2. We 
thus find a tower of resonances, corresponding to all positive values of n; the n-th resonance 
occurs at j n = n/2/i 2 , independently of K. If we further impose the Hamiltonian constraint, 
then K must take the value K n = (/ij n ) 2 = (n/2fi) 2 . 

It can be shown that the n/th resonance is surrounded by a stochastic layer whose width 
scales as n 1 / 3 /2/i 2 . For large n each resonance is much wider that the separation 1/2/z 2 
between resonances. Thus the stochastic layers merge and the behavior of the perturbed 
system becomes chaotic. 
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3.2. Chaos and Quantum Cosmology 

Though the ubiquity of chaos in cosmological models is interesting for its own sake, for the 
purposes of this review the point is how classical chaos impacts on quantum cosmological 
models. To analyze this issue, we shall show that a WKB solution for small Universes 
may diverge from the exact solution as the Universe expands, thus making the WKB 
interpretation of quantum mechanics unviable as is (61] EH I195[ 1233] . 

To see in simple terms why a WKB solution of the Wheeler - DeWitt equation departs 
from the exact solution when the underlying classical cosmological model is chaotic, let us 
consider again the procedure whereby the WKB solution is obtained in a multidimensional 
mini super space model. 

Let us return to a generic mini super space model as described in equations EH [29] and 
l30l Adopting the "covariant operator ordering, the Wheeler - DeWitt equation becomes 

H ' = {^^^W + ^} V = (54) 

where / = det f a p. We regard f as a scalar defined on the mini super space manyfold. We 
postulate a solution of the form 

= A 1/2 e im i s (55) 

We ask that [Ml holds order by order in mi. In the limit mi — > oo the dominant terms are 
O (jnp^j and 0(1). The first yields 

2 J dq a dqP m V ; 

and the second 

We understand [56] as meaning the more general equation 

^§W + ul ^ E (58) 

plus the constraint E = 0. We recognize [58] as the Hamilton- Jacobi equation for the mini 
super space evolution. The general solution will take the form S = S [q a , J a], where the J a 
are parameters regarded as world scalars. In general E = E [J a], only configurations such 
that E = E [J a] = are physical. 

Observe that [57] admits the solution [2211 269J 
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To reveal the geometric content of 
transformation to new variables 
dS 



e- 



and then 



dq 1 



observe that we may use S to generate a canonical 

(60) 

(61) 



The WKB wave function will blow up everytime the matrix inside the determinant becomes 
singular, that is, at caustics of the classical evolution. In a regular evolution these caustics 
are isolated; in a chaotic evolution caustics proliferate and the WKB wave function becomes 
arbitrarily "spiky". The exact solutions of the Wheeler - De Witt equation, on the other 
hand, are regular and cannot reproduce this behavior, therefore the WKB and exact solutions 
diverge from each other [321 1232] . 

This effect has been demonstrated through numerical solutions of the Wheeler - DeWitt 
equations and their WKB approximations for the conformal field FRW mini super space in 
refs. [611 and l89l. 



4. Decoherence and Cosmology 

4-1. Environment induced decoherence 

In this review, we shall stick to the position that the quantum mechanics of a closed system is 
unitary, and that therefore any decoherence effect must be caused through the interaction of 
the system proper with an environment. As we have remarked in the Introduction, while the 
formalism itself seems to place no requeriments on which systems qualify as environments, 
in practice not all choices are equally compelling. 

In the discussion of classical equations from quantum dynamics, Gell-Mann and Hartle 
[119J pointed out that for a large and possibly complex system the variables that will become 
classical habitually are the local densities of conserved hydrodynamic variables integrated 
over small volumes. To show that some variables become classical involves showing that they 
are readily decoherent and that they obey deterministic evolution equations. Hydrodynamic 
variables such as energy, momentum and number are of such characters because they are 
conserved quantities [U EJ E3 HU 11341 1135] . In [72] it is hown that these variables obey 
deterministic equations in spite of being strongly coupled to the environment because their 
inertia is correspondingly very high 

When such a natural choice is not forthcoming, one possibility is to choose an 
environment that warrants stability and predictability of the system proper dynamics. This 
leads to choose the environment by employing the Paz and Zurek's predictability sieve 
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[244, I288[ 1290] . The predictability sieve is the most general criterium available for the 
selection of the preferred basis. It has ben applied to cosmology in [l$3j 

For simplicity, we shall only consider cases where system proper and environment are 
actually two different systems. We shall only consider histories where the path ip a {t) of the 
system proper is specified, while the path of the environment is left unconstrained. 

If we assume an uncorrelated initial condition, then computing the decoherence functional 
reduces to a path integral over the environment 

D(a,a') = p^^A 8 '^- 8 -^^"^'^ (62) 

where S s is the free action of the system, p s is the initial density matrix for the system, and 
S IF is the Feynman - Vernon influence action [47J HH 11091 IllOj I124j 

jM^*) = f d& dc de, f 1 De f' De PE ^o s s ^- 8 ^^ ^ 

where Se is the action of environment, including the system-environment interaction and and 
P e is the environment initial density matrix. The weak decoherence condition is recovered 
when 

e-miSirWW)] « i Im[S IF {i)\i) 2 )}»\ (64) 

We observe that the influence action S/f = T is also the Schwinger - Keldysh (or closed 
time path - CTP) effective action for the environment under the approximation that system 
variables are treated as c-numbers, that is, Feynman graphs containing system propagators 
as internal lines are neglected [H [151 E3J M, ES E21 QZ21 Q2U ESS [2621 ESI]. This analogy 
implies that if we seek the effective dynamics for the system by looking at the stationary 
phase points of the decoherence functional, then this dynamics will be causal and dissipative. 
On the other hand, we may use the Feynman - Vernon trick of deriving the influence action 
from coupling the system to a random external source; then we see that the system dynamics 
is stochastic as well. 

The physical basis of dissipation and noise in the system dynamics is discussed in detail 
in [73]. What is going on is that the time evolution of the system variables excites the 
environment through parametric amplification. System and environment become entangled, 
and thus tracing over the second decoheres the first. The energy for the environment 
excitation comes from the system, whose own dynamics is therefore damped. The back 
reaction from the environment on the system has a random component (in the absense of 
special correlations in the initial state) which is perceived by the system as noise. The 
unitarity of the underlying theory enforces certain relationships between dissipation and 
noise, in an equilibrium situation these are just the fluctuation - dissipation theorem. 

We shall now see how the above applies to quantum cosmology. 
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4-2. Midi Super Space models 

The basic setup to discuss decoherence in quantum cosmology is to assume the mini 
superspaces we have studied above are simply the relevant sector of the theory, thus defining 
everything else in super space as the environment [ElO [1581 H521 ESQl MH UM IMS M31 [258] . 
However, since we do not really know how to do quantum cosmology (in the Wheeler - 
DeWitt sense) in super space, some further simplification is called for. A common approach 
is to reduce the number of degrees of freedom at the classical level, obtaining a theory with 
still infinite degrees of freedom but no general dipheomorphism invar iance, and thus more 
amenable to quantization. The resulting theory is called a midi super space. In this review 
we shall consider midi super spaces consisting of a Friedmann - Robertson - Walker Universe 
plus homogeneous scalar matter and linearized fluctuations around them. There are still 
many midi super spaces answering to this description, as we may consider all scalar, vector 
and tensorial perturbations, and then further different ways of splitting them into relevant 
variables and environment. For example, we may consider scalar fluctuations as relevant 
and tensorial as environmental, or the scalar fluctuation occupation numbers as relevant and 
the mode phases as environment, or some scalar fluctuations as relevant and the rest as 
environment, and so on [391 EH EHJ EQj MB EMl MM CEHEl EHE1 M M HHHl MM EUl MB 
[2071 12221 [2231 [259] . 

A basic discussion of decoherence, dissipation and noise in a FRW minisuperspace 
coupled to a massive conformally coupled free scalar field is given in [68]. However in 
this work linearized gravitational fluctuations were totally neglected, to the effect that 
the momenta constraint were violated. This may raise doubts about the soundness of the 
analysis. We shall now show that a more detailed analysis employing so-called gauge invariant 
variables [3 El EGI1 [1221 E25J [2261 [223 EEL] yields essentially the same results. If anything, 
the effects of noise and dissipation have been underestimated. 

For the sake of discussion we shall consider only the conformal factor of a spatially closed 
Friedmann - Robertson - Walker Universe and a homogeneous, conformally coupled massive 
real scalar field as relevant variables, and the scalar linearized fluctuations around them 
as the environment. In principle there are four different scalar fields associated to metric 
fluctuations, as we may perturb the lapse N = 1 + A, the shift A^ a = ^B^, the conformal 
factor a — )■ a + 5a = a (1 + (1/2) -0) or the three-metric g a b = 1 + E\ ah — (1/3) jabAE. There 
is one more scalar field associated to the matter fluctuations $ = + ip. Here j a b is the 
natural metric of a three sphere in the corresponding coordinates. We must also decompose 
the canonical momenta associated to the conformal factor P = Pq + 5P and matter field 
p = Po + Sp. The other geometrodynamic momenta are expressed in terms of a single scalar 
momentum conjugated to E. These make eight degrees of freedom per point. 

Neither all of these fluctuations are independent, nor they are physical. Upon a 
coordinate change x M — > x'^ = + £ M the metric changes into g'ffl = — — ^ i(U 
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and the field into $' = $ — £ M $ !At . The corrections to the background quantities are scalar 
perturbations provided £ a = £ >a . A crucial observation is that these "gauge" transformations 
mix the canonical variables with the lapse and shift variables |264] . Since we have two 
scalar functions £ and £o to play with, we may impose two arbitrary conditions. We choose 
B = E = 0. If we are given a coordinate system where these conditions do not hold, we may 
perform a coordinate change as above with £ = a 2 E/2 and £o = a 2 (B — Ej/2). This means 
that instead of the original fluctuations, we take as degrees of freedom the fluctuations in the 
frame where B = E = 0. Since this frame is (at least locally) unique, the resulting degrees 
of freedom are "gauge invariant" . 

The two "gauge fixing conditions" B — E — play formally the role of two new 
primary constraints on the theory, and so there also bring two secondary constraints meaning 
B = E = 0. Together with the two true secondary constraints H° = H a i a = we have in 
total six conditions the linearized fluctuations must satisfy, bringing the number of true 
degrees of freedom to just two. Moreover, the resulting theory may be brought through a 
canonical transformation into the theory of a canonical scalar field with a variable mass (see 
below). The difference with respect to the analysis in [68J is that the mass depends not only 
on the conformal factor (as is the case in the classical theory) but also on the homogeneous 
matter field. In other word, the time-dependent homogeneous matter field and the canonical 
gauge invariant variable are gravitationally coupled. Therefore, not just the expansion of 
the Universe, but matter field fluctuations too, may bring forth parametric amplification of 
the canonical gauge invariant field. 

Let us expand on the relationships among the linearized fluctuations. Observe that if 
B = E = then K a t = and then 7r a 5 = 0. This in turn implies E tt = 0, so the gauge 
condition is consistent. To check whether 71"^ = at the linearized level, we must compute 
the variation of the Hamiltonian w.r.t. variations 5 a b with j^dab = 0. This implies that 
g = 7 is not varied. 

Recall that on the background we have R^ cd = S^Jm — ^dlbc, Rab = ^lab and R = 6, 
Therefore SR = 2^ hd ^ ac 8(ic\ba ■ The variation of the Hamiltonian reads 



We must still consider the constraints of the theory. The scalar momentum constraint reads 




(65) 




(66) 



(Fo, a +p$,,)-^ ^(Pa+p®) 







(67) 
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Linearizing we get 

P aip + 2p {p — a5P — (f)5p = 
and the Hamiltonian constraint 



(68) 



1 



p2 + 



--P 



24m2^7 2^7 
2m 2 p ^{-Aa 2 + 3r\aa b } 



+ 2^ 



(69) 



Whose linear part reads 



1 1 
——P 5P + -p 5p - 6m 2 „a 2 ij + 0y? + 2m 2 a 2 Aip 



+ 



rrr 



2 

'■a 2 (<p 2 i> + 20v?) - -<j)Aip 







(70) 



We may use both constraints to express ip and 5P as functionals of ip and 5 p , or viceversa. 
If we make the first choice, then the dynamics of ip and 5p is described by an effective 
Hamiltonian density 



H, 



eff 



H - SP5a 



We read Sa tt out of the original Hamilton equations 

1 



5a, t = 



[5P + AP 



(71) 



(72) 



We also expand the scalar linear fluctuations into eigenfunctions of the spatial Laplace 
operator 

J2 h (*) (t);5p=^J2 h (x) 5 PL (t) (73) 



Af L = —L 2 fi 



(74) 



where L represents the integers L, K, M with L > 0, —L < K, M < L, and L 2 = L (L + 1). 
We assume the Jl are orthonormal. Then the linearized equations for different modes 
decouple. Solving the linearized equations for each mode, all perturbations may be expressed 
in terms of <£l and 8pi- The effective Hamiltonian becomes a sum over modes, each mode 
being a squeezed harmonic oscillator. These may be reduced by standard methods to 
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the Hamiltonian of a harmonic oscillator with a time-dependent instantaneous frequency 
[73] . The frequency will depend on the scale factor, the homogeneous scalar field and their 
derivatives, though the latter dependences shall be suppressed by powers of the Planck mass. 
Therefore not only the expansion of the Universe, but also a time-dependent matter field, 
will lead to parametric amplification of the fluctuations. 

4-3. Decoherence and particle creation 

We shall now carry out an analysis of noise, fluctuations and dissipation in the 
minisuperspace as induced by the gauge invariant fluctuations. 

As we have seen, each mode of the gauge invariant scalar variable may be reduced to a 
parametric oscillator. Quantization proceeds by further decomposing each Fourier amplitude 
in its positive and negative frequency parts, defined by a suitable choice of time parameter. 
This is accomplished by developing the corresponding mode on a basis of solutions of the 
Klein - Gordon equation, so normalized that the positive frequency function has unit Klein 
- Gordon norm, the negative frequency function has norm —1, and they are mutually 
orthogonal in the Klein - Gordon inner product. Such a basis of solutions constitute a particle 
model. Properly normalized particle models are related to each other through Bogolubov 
transformations. Further identification of the coefficient of the positive frequency function 
in the development of the field, as a destruction operator, allows for the second quantization 
of the theory. The particle model is also associated to a vacuum state, which is the single 
common null eigenvector of the destruction operators, and to a Fock basis, built from the 
vacuum through the action of the creation operators [31 [2271 HH 12351 12351 1237J 12351 1239] 

It is also well known that in a generic dynamic space time, there is no single particle 
model which can be identified outright with the physical concept of "particle"; however, 
oftentimes it is possible to employ a variety of criteria (such as minimization of the particle 
number as detected by a free falling particle detector, Hamiltonian diagonalization, conformal 
invariance, analytical properties in the euclidean section of the space time, if any, etc) to 
single out a preferred particle model in the distant past (or "IN" particle model) [56] [139]. 
and another in the far future, or "OUT" particle model. In general, these models are not 
equivalent, the vacuum of one model being a multiparticle state in the other. In our problem, 
the choice of boundary conditions for the path integral above amounts to a definite choice 
of the IN particle model, and the IN quantum state, in each branch of the closed time path 

m 

Since the CTP effective action is independent of the OUT quantum states, we have 
more freedom in choosing an OUT particle model. It is convenient to choose a common 
OUT particle model for both evolutions (that is, the Cauchy data on the matching surface 
1] = r]° are the same although the actual basis functions will be different). The positive- 
frequency time dependent amplitude functions f± for the conformal model in each branch 
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are related to those F of the OUT model by f± = a±F+/3±F* at 77 = 77 , where oc±, (3± are the 
Bogolubov coefficients in each branch, obeying the normalization condition |a±| 2 — |/3±| 2 = 1. 
The CTP effective action is found to be [HHl E3] : 

r = (|)ln[a_<-/3_/3;] (75) 

This expression is an exact evaluation of the Gaussian path integral defining the influence 
action. 

If we regard the saddle points of the decoherence functional as the actual classical 
trajectories, we see that the decoherence functional is essentially the Feynman - Vernon 
influence action for the minisuperspace variables interacting with the environment provided 
by the gauge invariant fluctuations; the influence functional is also closely related to the 
Schwinger - Keldish effective action. By following through with this analogy, it is easy to 
show that the resulting equations of motion are real and causal. 

The lesson for us is that there can be decoherence (ImT > 0) if and only if there 
is particle creation in different amounts in each evolution. (This is also implicit in 
[242[ 1243] . ) Indeed, we can always choose the OUT model so that a + = 1, + = 0, yielding 
T = (i/2)lno!_. The condition for decoherence in this case is then |a_| > 1. But since 
I or_ 1 2 = 1 + | 2 , this can only happen if there is particle creation between these two 
particle models. 

To quantify the above considerations, let us return to the FRW plus a massive conformal 
field mini super space. Let us consider the system evolution to be close to the n-th 
resonance. This means Uj = n/2, uk = —1/2. We assume n ^> 1 and m < m p . 
Since the field oscillations are much faster than the expansion and contraction of the 
Universe we may freeze the conformal factor at its peak value. Therefore in the background 
a = n/m and Pq = 0, while, with a convenient choice of phase, = vl2 (m p /m) sin nrj 
and p = 2nix 2 \/Y2 (m p /m) cosnrj. Under these conditions the gauge invariant fluctuations 
become A = ip = and SP = (2p /a) (p = An 2 \f\2m p (cos nrj) (p. Therefore the L-th mode 
<Pl is a periodically driven oscillator with instantaneous frequency fl 2 = + ufcos2nr/ [73] 
with u 2 = n 2 + L 2 + L\ [L\ = 1 + 8tt 4 ) and u\ = 8tt 4 . Because ojq 3> u\ this is a very narrow 
resonance and the effects of parametric amplification will not be strong. 

However, when we include chaos in our considerations the overall effect is that the 
instantaneous frequency for each mode acquires a random component. This situation has 
been studied by Zanchin et al |278[ 1279] . who show that randomness amplifies particle 
creation in all modes. Therefore we are led to the same conclusion than |198j . namely, the 
relevant parameter range for chaos is also the most favourable one for particle creation, and 
therefore for decoherence. 

The physical mechanism underlying decoherence in the decoherent history scheme of 
Gell-Mann and Hartle is the same as in the environment-induced scheme based on a reduced 
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density matrix obtained by tracing over the environmental degrees of freedom, namely the 
entanglement of system and environment through particle creation |^9 | 150 ] [57 ] E5 j I161[ f!62j. 
Since the energy of the created particles must be supplied by the minisuperspace variables, 
it is clear that the back reaction of the environment on the system proper will have a 
damping effect. Moreover, quantum fluctuations in the number of created particles result in 
fluctuations in the back reaction, which are perceived by the system proper as a stochastic 
driving [73J. The joint effect of dissipation and noise may induce a qualitative change in 
the dynamics of the system proper, which is most conspicuous when this system is chaotic 
[198J. Cosmological models where this stochastic driving is explicitly taken into account are 
discussed in [Ml Ml EE1 ED ESI USB ED UM EDH1 EDI 1205] . 

4-4- Chaos, Decoherence and Quantum Cosmology 

In ordinary quantum mechanical systems, the best way to visualize the effect of noise and 
dissipation on quantum chaos is the Wigner function |277] . Recall the the dynamics of a 
chaotic system, for example in the neighborhood of an equilibrium state, is exponentially 
contracting in some directions in phase space, while expanding along others. Up to terms 
of order h 2 the Wigner function obeys the classical Liouville transport equation. Therefore 
when the system is chaotic, it tends to produce structure in arbitrary short scales along the 
contracting directions of the dynamics, while smoothing out along the expanding directions. 
However, when noise and dissipation are factored in, the resulting diffussion introduces a 
shortest scale where structure may be produced. On the other hand, the expansion of the 
Wigner function along the unstable directions is limited by the Liapunov exponents, closely 
related to the Kolmogorov entropy through Pesin's theorem. The combination of these two 
effects implies that the volume of the set where the Wigner function is substantially non zero 
increases in time, until it saturates. Since the logarithm of this volume is the thermodynamic 
entropy, we find that the thermodynamic entropy of the system proper is increasing at a 
rate which is determined by the Kolmogorov entropy of the system itself. Observe that 
the rate of entropy increase is largely independent of the system - environment interaction. 
This simple but compelling picture of the effect of decoherence on quantum chaos has been 
put forward by Paz and Zurek [291] and since confirmed by detailed numerical simulations 

[21 El HH EH1 EH rrrni H2H1 H2HI EIHi [2201 E57]. 

To see whether this scenario obtains in quantum cosmology recall that in quantum 
cosmology we regard the Wave function of the Universe as a scalar field on minisuperspace. 
The density matrix is therefore akin to a Schwinger function, namely the expectation value 
of the product of the field evaluated at two different locations. If the minisuperspace is flat, 
then the Wigner function may be introduced as a Fourier transform of the density matrix 
with respect to the difference between its arguments. In the general case, however, mini 
super spaces are described by curved manifolds. In this general case, notwithstanding, a 
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local momentum space may be introduced by using the same techniques as applied to derive 
quantum kinetic theory from quantum field theory on curved spaces [62 | [65 | [73 |ll26[ll27j [128j. 

The Wigner function will generally obey two equations, a "mass shell condition" and 
a first order dynamic equation. Both equations are covariant in phase space. In the WKB 
limit, a well defined concept of time arises, namely the time of the classical solution around 
which the fluctuations are defined. In this limit the dynamic equation becomes equivalent 
to the usual equations investigated by Paz and others. 

Therefore it seems appropriate to conclude that noise and dissipation due to parametric 
excitation of gauge invariant fluctuations restores the validity of the WKB limit by avoiding 
the proliferation of spikes in the Wave Function of the Universe, and produce both 
decoherence and entropy production at a rate determined largely by the mini super space 
dynamics. However, it is fair to say that this issue has not been explored in the literature with 
the detail it demands, and therefore our arguments do not intend to be taken as anything 
but strong hints 

5. Final Remarks 

In this review we have briefly recalled forty something years of efforts towrads the unraveling 
of the relationships between chaos, decoherence and quantum cosmology. The result 
notwithstanding seems patchy and inconclused, because some obvious leads have not been 
thoroughly followed up yet. The most glaring gap, in my view, is whether the translation 
of the sustantial body of work on Wigner functions for classically chaotic systems may hit 
some snag due to the specifics of reparametrization invariant theories. To the best of my 
knowledge this issue has not been investigated beyond the semiclassical approximation. 

Given the size of the reference section, it may be useful to conclude by pointing out 
a few highlights; this is of course a very subjective choice on top of the arbitrary selection 
which created the reference list in the first place. 

To begin with, any discussion of quantum cosmology should be set within the framework 
of our current understanding of the Very Early Universe. Mukhanov's |225j stands out among 
many excellent books and is also a very good introduction to the formalism of gauge invariant 
variables. 

Besides the foundational work such as [92j I146j 1272] and [276], references [133] and 
[182J are excellent entry points to all aspects of quantum cosmology. Similarly, the books 
by Giulini et a/ |121] . Schlosshauer [253] and Weiss [275] are well known introductions to the 
physics of quantum open systems. The classic exposition of the Gell-Mann/ Hartle formalism 
is [TT7] . 

The foundational papers of cosmological chaos [26] and [213] are undisputed 
masterpieces. Together with the review [20] they cover much of what is there to be said. For 
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an updated state of the art see |152j . 

With regards of chaos in FRW cosmologies, the early paper [58] already raised many 
relevant issues. For the state of the art see (HSJ I201[ [202J. 

The 1993 Kananaskis meeting |156j was arguably a turning point in the development 
of the subject. Although I have quoted several individual papers from its Proceedings, the 
volume clS 8u whole has a remarkable sinergia. 

To the best of my knowledge, (61] was the first detailed demonstration of the divergence 
between solutions of the Wheeler - DeWitt equation in a mini super space and their WKB 
approximations. See also [89] . 

Our discussion of decoherence in quantum cosmology was essentially the application to 
the quantum cosmological context of tools which are well proven in the much wider fields 
of environment induced decoherence and nonequilibrium quantum field theory. For more 
general introductions to these subjects see [245] and [73] . respectively. References [167J and 
[68] may be regarded as blueprints for a decoherent quantum cosmology program. For the 
further progress of the stochastic cosmology program see [171] . 

In spite of the strong theoretical foundations [32] and ground breaking work by 
Habib and others [1261 11271 1128] , the application of Wigner functions to classically chaotic 
cosmological models is still underdeveloped. I call on my readers to finish up the task. 
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